Time-Delayed Feedback in Neurosystems 
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The influence of time delay in systems of two coupled excitable neurons is studied in the framework 
of the FitzHugh-Nagumo model. Time-delay can occur in the coupling between neurons or in 
a self-feedback loop. The stochastic synchronization of instantaneously coupled neurons under the 
influence of white noise can be deliberately controlled by local time-delayed feedback. By appropriate 
choice of the delay time synchronization can be either enhanced or suppressed. In delay-coupled 
neurons, antiphase oscillations can be induced for sufficiently large delay and coupling strength. 
The additional application of time-delayed self-feedback leads to complex scenarios of synchronized 
in-phase or antiphase oscillations, bursting patterns, or amplitude death. 

1. INTRODUCTION 

Control of unstable or irregular states of nonlinear dynamic systems is a central issue of current research 
A particularly simple and efficient control scheme is time-delayed feedback which occurs naturally in a number of 
biological systems including neural networks where both propagation delays of electrical signals connecting different 
neurons and local neurovascular couplings lead to time delays [2, Moreover, time-delayed feedback loops might 

be deliberately implemented to control neural disturbances, e.g., to suppress undesired synchrony of firing neurons in 
Parkinson's disease or epilepsy 0, 0, 0]- Chaos control techniques have been first applied experimentally in vitro in a 
spontaneously bursting neural network Q using a method proposed by Ott, Grebogi and Yorke [8]. In contrast to this 
method, the time-delayed feedback method by Pyragas [§] and its extensions [TO] do not require detailed information 
on the system to be controlled, and they are non-invasive, i.e., once the target state is reached, the control force 
vanishes. However, one needs to tune the feedback gain and the delay time. Therefore, theoretical investigations help 
to understand optimal choices of these parameters. Various global delayed feedback schemes have been proposed as 
effective and robust therapy of neurological diseases with pathological synchronization causing tremor [y, |7| . They 
have been contrasted with local delayed feedback methods [Til . [T^ . \\^ . 

Neurons are excitable units which can emit spikes or bursts of electrical signals, i.e., the system rests in a stable 
steady state, but after it is excited beyond a threshold, it emits a pulse. In the following, we consider electrically 
coupled neurons. Such electrical synapses are less common than chemical synapses, but it has been shown in hip - 
pocampal slices that high-frequency synchronized oscillations are independent of chemical synaptic transmission [15j , 
and physiological, pharmacological, and structural evidence was provided [T^ that axons of hippocampal pyramidal 
cells are electrically coupled (Fig. [T] (a)). Time delays in the coupling must be considered particularly in the case of 
high-frequency oscillations. 

The simplest model that displays features of neural interaction consists of two coupled neural systems. A single 
electrical synapse can lead to synchronized co-operative behaviour between two axo-axonally coupled hippocampal 
pyramidal cells when only one of them is stimulated antidromically at high frequency [la, [13 ■ Moreover, such 
neurons have been shown to reflect antidromic spikes (propagating opposite to the normal direction) at the soma, but 
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FIG. 1: (a) Scheme of two axo-axonally coupled neurons (pyramidal cells coupled by an electrical synapse) [3l. (b) Two 
mutually coupled neural populations (delay t, coupling constant C) with feedback control loop (delay tk, coupling constant 
K) and noise input Di, D2- 
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a continuous reverberating activity is not generated by axonal reflections in two axo-axonally coupled cells alone [lj| . 
In order to describe the complicated interaction between billions of neurons in large neural networks, the neurons are 
often lumped into highly connected sub-networks or synchronized sub-ensembles. Such neural populations are usually 
localized spatially and contain both excitatory and inhibitory neurons ^1^. In this sense, the model of two mutually 
coupled neurons may also serve as a paradigm of two coupled neural sub-ensembles. 

Starting from such simplest network motifs, larger networks can be built, and their effects may be studied. For 
example, starting from two interconnected reticular thalamic neurons with oscillatory behavior, it was shown in [Toj 
how more complex dynamics emerges in ring networks with nearest neighbors and fully reciprocal connectivity, or in 
networks organized in a two-dimensional array with proximal connectivity and "dense proximal" coupling in which 
every neuron connects to all other neurons within some radius. In another example [20|, a neural population was 
itself modeled as a small sub-network of excitable elements, to study hierarchically clustered organization of excitable 
elements in a network of networks. 

In the following, we consider two mutually coupled neurons modelled by the paradigmatic FitzHugh-Nagumo system 
pll 113, H^l in the excitable regime: 

2-3 

eiXi = xi - ^ - yi + C[x2{t - t) - xi{t)] 
yi = xi +a + Di^i(t) 

62X2 = X2 - ^ - y2 + C[xi{t - t) - X2{t)] 

2/2 = X2+a + D2^2it) (1) 

where xi, yi and X2, 2/2 correspond to single neurons (or neuron populations), which are linearly coupled with coupling 
strength C. The variables xi, X2 are related to the transmembrane voltage and yi, y2 refer to various quantities 
connected to the electrical conductance of the relevant ion currents. Here a is an excitability parameter whose value 
defines whether the system is excitable (a > 1) or exhibits self-sustained periodic firing (a < 1), ei and 62 are the 
timescale parameters that are usually chosen to be much smaller than unity, corresponding to fast activator variables 
Xi, X2, and slow inhibitor variables j/i, 2/2- 

The synaptic coupling between two neurons is modelled as a diffusive coupling considered for simplicity to be 
symmetric [IJ, [Ig] . More general delayed couplings are considered in [23| • The coupling strength C summarizes 
how information is distributed between neurons. The mutual delay r in the coupling is motivated by the propagation 
delay of action potentials between the two neurons xi and X2- 

Each neuron is driven by Gaussian white noise ^i{t) {i — 1,2) with zero mean and unity variance. The noise 
intensities are denoted by parameters Di and D2, respectively. 

Besides the delayed coupling we will also consider delayed self-feedback in the form suggested by Pyragas , where 
the difference s{t) — s(t — t^) of a system variable s (e.g., activator or inhibitor) at time t and at a delayed time 
t ~ tk, multiplied by some control amplitude K , is coupled back into the same system (Fig. (Ub)). Such feedback 
loops might arise naturally in neural systems, e.g., due to neurovascular couplings that has a characteristic latency, 
or due to finite propagation speed along cyclic connections within a neuron sub-population, or they could be realized 
by external feedback loops as part of a therapeutical measure, as proposed in Ref. Q. This feedback scheme is 
simple to implement, quite robust, and has already been applied successfully in a real experiment with time-delayed 
neurofeedback from real-time magnetoencephalography (MEG) signals to humans via visual stimulation in order to 
suppress the alpha rhythm, which is observed due to strongly synchronized neural populations in the visual cortex in 
the brain [281] . One distinct advantage of this method is its non- invasiveness, i.e., in the ideal deterministic limit the 
control force vanishes on the target orbit, which may be a steady state or a periodic oscillation of period r. In case 
of noisy dynamics the control force, of course, does not vanish but still remains small, compared to other common 
control techniques usingexternal periodic signals, for instance, in deep-brain stimulation to suppress neural synchrony 
in Parkinson's disease }29| . 

The phase portrait and the null-isoclines of a single FitzHugh-Nagumo system without noise and feedback are 
shown in Fig. [2l[a). The fixed point A is a stable focus or node for a > 1 [excitable regime). If the system is perturbed 
well beyond point A' (see inset), it performs a large excursion A^B^C^D^A'va phase space corresponding 
to the emission of a spike (Fig. Hfb)). At a = 1 the system exhibits a Hopf bifurcation of a limit cycle, and the fixed 
point A becomes an unstable focus for a < 1 [oscillatory regime). 

In the following we choose the excitability parameter a = 1.05 in the excitable regime close to threshold. If noise 
is present, it will occasionally kick the system beyond A' resulting in noise-induced oscillations [spiking). 
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2. STOCHASTIC SYNCHRONIZATION OF INSTANTANEOUSLY COUPLED NEURONS 

We shall first consider two coupled FitzHugh-Nagumo systems as in eq. jTOj albeit without delay in the coupling 
(r = 0). Noise can induce oscillations even though the fixed point is stable [30l[3l|. The noise sources then play the 
role of stimulating the excitable subsystems. Even if only one subsystem is driven by noise, it induces oscillations of 
the whole system through the coupling. In the following, we consider two nonidentical neurons, described by different 
timescales ei = 0.005, £2 = 0.1, and set the noise intensity D2 in the second subsystem equal to a small value, 
D2 = 0.09, in order to model some background noise level. Depending on the coupling strength C and the noise 
intensity Di in the first subsystem, the two neurons show cooperative dynamics. 

Figure [3l^a) depicts the temporal dynamics of a single FitzHugh-Nagumo system (C = 0) due to stochastic input. 
One can see that the system is excitable since it performs large excursions in phase space. Figure Eljb) shows the 
temporal dynamics of the two instantaneously coupled neural systems for increasing noise intensities Z?i, where the 
green (light grey), red (dark grey), and black curves correspond to xi, X2, and their sum = xi + X2, respectively. 
For Di = 0, the first subsystem is enslaved and emits a spike every time the second unit does. For increasing Z3i, 
however, this synchronization is weakened and for large noise intensities Di the dynamics of the first subsystem is 
independently dominated by its own stochastic input. 

There are various measures of the synchronization of coupled systems j33 [. For instance, one can consider the 
average interspike intervals of each subsystem, i.e. (Ti) and {T2), calculated from the x variable of the respective 
subsystem. Their ratio {Ti)/{T2) is a measure of frequency synchronization, as depicted in Fig. IH Panel (a) displays 
the ratio of average interspike intervals in dependence on the noise intensity Di (green dots) for fixed coupling 
strength C = 0.07. One can see that for increasing Di the ratio {Ti) / {T2) decreases. Thus, the two subsystems 
become less synchronized. Panel (b) shows the dependence on the coupling strength C for fixed noise intensity 
Di = 0.25. Without coupling (C = 0), the two subsystems operate on their own timescale. For increasing C, however, 
they synchronize as the ratio of the average interspike intervals approaches unity. Panel (c) shows the result as a 
function of Di and C, where the bright regions indicate a strongly synchronized behaviour of the two subsystems. 
For small Di and large coupling strength C, the two subsystems display synchronized behaviour, {Ti) / {T2) ~ 1. On 
average, they show the same number of spikes indicated by bright (yellow) color. In the next Section, we consider 
three different regimes of synchronization: moderately, weakly, and strongly synchronized systems as marked by 
black dots in Fig. [4l[c). 

Other measures for stochastic synchronization are given by the phase synchronization index ^30l | , or the mean phase 
synchronization intervals [33| . but they exhibit qualitatively similar behaviour. 
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FIG. 2: Excitable dynamics of a single FitzHugh-Nagumo system: (a) Phase portrait (2:1, 3/1) (trajectory: solid blue, nuUclines: 
dashed black), (b) Time series of activator xi{t) (red) and inhibitor yi{t) (green). The colored dots A, B, C, D mark corre- 
sponding points on panels (a) and (b). The inset in (b) shows a blow-up of the phase portrait near A. Parameters: ei — 0.01, 
a = 1.05, Di = 0. 
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FIG. 3; (a): Noise-induced oscillations in a single FitzHugh-Nagumo system (C = 0) for a noise intensity Di = 0.6. The red 
dashed and dash-dotted lines denote the null- isoclines. Panel (b) : Time series of two instantaneously coupled FitzHugh-Nagumo 
systems for different noise intensities Di. The green (light grey), red (dark grey), and black curves correspond to xi, X2, and 
their sum xs ~ xi + X2, respectively. Other parameters: a = 1.05, ei = 0.005, €2 = 0.1, r = 0, C = 0.07 and D2 — 0.09. 




FIG. 4: Frequency synchronization: Average interspike intervals of each subsystem (black solid curve: (Ti), red dashed: (T2)) 
and their ratio (green dotted: {Ti) / {T2)) (a) vs. noise intensity Di (for fixed C = 0.07), (b) vs. coupling strength C (for fixed 
D\ = 0.25). (c): Ratio of the average interspike intervals in the (Di, C) plane. The black dots mark the parameter choices for 
different synchronization regimes used in Fig. [S] Other parameters: a — 1.05, ei — 0.005, £2 ~ 0.1, t — 0, and D2 = 0.09. 



3. CONTROL OF SYNCHRONIZATION BY TIME-DELAYED FEEDBACK 

In this Section, we consider the control of global cooperative dynamics by local application of a stimulus to a single 
system. This stimulus is realized by time-delayed feedback, which was initially introduced by Pyragas in order to 
stabilize periodic orbits in deterministic systems Q: 

eiii = a;i - y - j/i + C[a;2(<) - xi(t)] (2a) 
yi = xi + a + K [yi{t - tk) - yi{t)] + Di^i{t) 

£2x2 = a;2-y-?y2 + C[a;i(t)-X2(t)] (2b) 

2)2 = X2+a + D2^2{t). 
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The parameters of the time-delayed feedback scheme are the feedback gain K and the time delay tk- With this 
method, a control force is constructed from the differences of the states of the system which are one time unit 
apart. One could also consider application of the feedback scheme to both subsystems and effects of different values of 
the control parameters for each subsystem, but these investigations are beyond of the scope of this work. Previously, 
time-delayed feedback has also been used to influence noise- induced oscillations of a sing le excitable system 34, 35, SGj], 
of systems below a Hopf bifurcation [53, [H, [s^ S^l or below a global bifurcation [4lL fi^ , and of spatially extended 
reaction-diffusion systems [l^, 13, S] • Extensions to multiple time-delay control schemes have also been considered 
[31, 33, 45, 4a] . 




FIG. 5: Control of synchronization: Interspike intervals vs. time delay tk- Green dots: ratio of the interspike intervals 
{T\) / {T2) of the two subsystems, black solid curve: (Ti), red dashed curve: (T2). Panels (a), (b), and (c) correspond to the case 
of moderate (C = 0.2, D\ = 0.6), weak (C = 0.1, D\ = 0.6), and strong (C = 0.2, Di — 0.15) synchronization, respectively. 
Other parameters: a = 1.05, ei = 0.005, £2 = 0.1, C = 0.07, r = 0, D2 = 0.09, K = 1.5. 




FIG. 6: Ratio of interspike intervals {T\) / {T2) in dependence on the feedback gain K and the time delay r for (a) moderate 
[C = 0.2, Di — 0.6), (b) weak (C = 0.1, Di = 0.6), and (c) strong (C — 0.2, Di = 0.15) synchronization. Other parameters as 
in Fig. SJc). 



Since we are interested in the effects of a control force on the synchronization, in the following we consider three 
different cases: moderately, weakly, and strongly synchronized systems given by the specific choices of the coupling 
strength and noise intensity in the first subsystem. These different cases of stochastic synchronization are marked as 
black dots in Fig. [H 

As a measure to quantify changes in the synchronization due to the control force, we consider the ratio of average 
interspike intervals. In the presence of a control force, i.e., K ^ 0, the cooperativity can be influenced by varying the 
feedback gain K and the time delay tk- 

For fixed feedback gain K = 1.5, Fig. [5] depicts the average interspike intervals of the subsystems, shown as solid 
(black) and dashed (red) curves for (Ti) and (T2), and their ratio (green dots) for the case of (a) moderately, (b) 
weakly, and (c) strongly synchronized systems, respectively, in dependence on the time delay tk- In all three cases, 
the stochastic synchronization can be strongly modulated by changing the delay time, i.e., one can either enhance 
and suppress synchronization by appropriate choice of the local feedback delay. 
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The overall dependence of the frequency synchronization, measured by the ratio of (Ti)/(T2), is displayed in 
dependence on the control parameters K and tx in Fig. [S] for moderate, weak, strong synchronization in panels (a) , 
(b), and (c), respectively. Thus, Fig.[5]can be understood as a horizontal cut for K = 1.5 through Fig.[6l One can see a 
modulation of the ratio of average interspike interval by tk for a large range of feedback gain. In view of applications, 
where neural synchronization is often pathological, e.g. in Parkinson's disease or epilepsy, it is interesting to note 
that there are cases where a proper choice of the local feedback control parameters leads to desynchronization of the 
coupled system (dark regions in Fig.[S]). 



4. DELAY-COUPLED NEURONS 



In this section we study the influence of a delay in the coupling of two neurons, rather than a delayed self- feedback. 
We set the noise terms in Eq. ([T]) equal to zero, Di = D2 = 0, but consider a time-delay r in the coupling. In the 
deterministic system the delayed coupling plays the role of a stimulus which can induce self-sustained oscillations in 
the coupled system even if the fixed point is stable. In this sense the delayed coupling has a similar effect as the 
noise term in the previous sections. Here the bifurcation parameters for delay-induced bifurcations are the coupling 
parameters C and r. 



4.1. Linear stability of fixed point 



In the following we shall choose symmetric timescales ti — €2 ~ e ~ 0.01 and fix a = 1.05, where each of the two 
subsystems has a stable fixed point and exhibits excitability. 

The unique fixed point of the system is symmetric and is given by x* = (x*, yj, X2, 2/2)1 where x* — —a, y* — /3~a. 
Linearizing Eq. ((T|) around the fixed point x* by setting x(i) = x* -I- (5x(<), one obtains: 
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where ^ = 1 



C. The ansatz 



(5x(t) 



where u is an eigenvector of the Jacobian matrix, leads to the characteristic equation for the eigenvalues A: 



(1 - + eX' f - {XCe-^^'f = 0, 



which can be factorized giving 



1 - ^A + eA^ ± ACe 



-At 



0. 



(4) 
(5) 
(6) 



This transcendental equation has infinitely many complex solutions A. Fig. [7] shows the real parts of A for various 
values of C. As can be seen in Fig. [7] the real parts of all eigenvalues are negative throughout, i.e., the fixed point 
of the coupled system remains stable for all C. This can be shown analytically for a > 1 by demonstrating that no 
delay-induced Hopf bifurcation can occur. Substituting the ansatz X = iuj into Eq. ^ and separating into real and 
imaginary parts yields for the imaginary part 



This equation has no solution for a > 1 since |^| 
occur. 



C = ±Ccos(wt) (7) 
= — 1 -|- C > C, which proves that a Hopf bifurcation cannot 



4.2. Delay-induced antiphase oscillations 



Delay-induced oscillations in excitable systems are inherently different from noise-induced oscillations. The noise 
term continuously kicks the subsystems out of their respective rest states, and thus induces sustained oscillations. 
Instantaneous coupling without delay then produces synchronization effects between the individual oscillators [sol [sil . 
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|33|. For delayed coupling the case is entirely different. Here the impulse of one neuron triggers the other neuron to 
emit a spike, which in turn, after some delay, triggers the first neuron to emit a spike. Hence self-sustained periodic 
oscillations can be induced without the presence of noise (Fig.[8|). It is evident that the oscillations of the two neurons 
have a phase lag of tt. The period of the oscillations is given by T = 2(r + 6) with a small quantity 5 > 0. 

In order to understand this additional phase-shift 6, we shall now consider in detail the different stages of the 
oscillation as marked in Fig. [21 Due to the small value of e <C 1 there is a distinct timescale separation between the 
fast activators and the slow inhibitors, and a single FitzHugh-Nagumo system performs a fast horizontal transition 
A B, then travels slowly approximately along the right stable branch of the xi nuUcline B C (firing), then 
jumps back fast to D, and returns slowly to the rest state A approximately along the left stable branch of the xi 

nullcline (refractory phase). If a is close to unity, these four points are approximately given hy A = (—a, —a + -^), 

B = (2, -|), C = (1, |), D = (-2, |). A rough estimate for A' is (a -2,~a + ^). The two slow phases B —> C and 

D ^ A can be approximated by yi w xi — ^ and hence j/i w ii(l — xf) — xi + a which gives 



Xi 



xi + a 



1 - xj 

which can be solved analytically, describing the firing phase (+) and the refractory phase (-): 



(8) 



dxi = (a^ - l)ln a(±2 - x) + 2 = t. 

j-2 xi + a a + X 2 



(9) 



Integrating from B to C gives the firing time 

r-l 

Tf 



dxi = (a'' - 1) In 

2 xi + a 



a + -. 

a + 1 2 



(10) 



For e = 0.01, a = 1.05 the analytical solution is in good agreement with the numerical solution in Fig.[2l^b), including 
the firing time Tf — 0.482 (analytical approximation: 0.491). 

For a rough estimate, in the following we shall approximate the spike by a rectangular pulse 



xiit) 



2 \i t<Tf 
-a a t> Tf 



(11) 



If the first subsystem is in the rest state, and a spike of the second subsystem arrives at t = (after the propagation 
delay r), we can approximate the initial dynamic response by linearizing xi,yi around the fixed point (x*,?/!!") and 



C = 0.1 



C = 0.4 




FIG. 7: Real parts Re(A) of the eigenvalues of the fixed point vs. time delay r for a = 1.05, e = 0.01, and (a) C=0.1, (b) C=0.4, 
(c) C=l, (d) C=2. 
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FIG. 8: Delay-induced oscillations, (a), (b): Time series of both subsystems (red solid lines: activator Xi, green solid lines: 
inhibitor yi\ black dashed lines: fixed point values of activator and inhibitor), (c), (d): Phase portraits of activators (c) and 
inhibitors (d). Parameters: a = 1.05, e = 0.01, C = 0.5, r = 3. 



approximating the feedback by a constant impulse during the firing time Tf. 
direction is then given by 

6(5x1 = i,6xi + 2C 



The fast dynamic response along the x\ 



(12) 



with ^ < 0. This inhomogeneous linear differential equation can be solved with the initial condition a;i(0) = —a: 

2C 



(1 



(13) 



Note that this equation is not valid for large t since (i) the linearization breaks down, and (ii) the pulse duration 
is exceeded. For small t Eq. p^ can be expanded as 



xi(t) 



2C 



-t 



(14) 



which is equivalent to neglecting the upstream flow field — |CI<^a;i in Eg. p^ near the stable fixed point A compared to 
the pulling force 2C of the remote spike which tries to excite the system towards B. Once the system has crossed the 
middle branch of the x\ nuUisocline at A', the intrinsic flow fleld accelerates the trajectory fast towards B, initiating 
the firing state. Therefore there is a turn-on delay 5, given by the time the trajectory takes from A to A', i.e. 
x\[6) a — 2, according to Eq. lM]) : 



(5 = (a - 1 



(15) 



Since the finite rise time of the impulse has been neglected in our estimate, the exact solution b is slightly larger and 
does not vanish at a = 1. 

With increasing a the distance A ~ A' increases, and so does S. The small additional phase shift S between the 
spike xi{t) and the delayed pulse X2{t — t) results in a non- vanishing coupling term at the beginning and at the end 
of the spike xi(t). It is the reason (i) that the spike is initiated, and (ii) that it is terminated slightly before the 
turning point of the xi nuUcline. The latter effect becomes more pronounced if a is increased or t is decreased (Fig. 
[HI) . Both lead to a shift of the initial starting point of the spike emission on the left branch of the nuUcline towards 
D, and hence to a longer distance up to the middle branch of the nuUcline which has to be overcome by the impulse 
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X2, hence to a larger turn-on delay 5, and therefore to an earlier termination of the spike Xi. This explains that the 
firing phase is shortened, and the limit cycle loop is narrowed from both sides with increasing a or decreasing r, see 
Fig. [9l In the case of a = 1.05 and t — i (Fig. [9] (a)), the delay time is large enough for the two subsystems to 
nearly approach the fixed point A before being perturbed again by the remote signal. If the delay time becomes much 
smaller, e.g., for r = 0.8 (Fig. [9] (b)), the excitatory spike of the other subsystem arrives while the first system is still 
in the refractory phase, so that it cannot complete the return D ^ A to the fixed point. In this case, a in Eq. p5|) 
has to be substituted by a larger value o with a < a < 1.7 in order to get a better estimate of 5. Note that without 
the phase-shift 5 the coupling term C[x2{t — t) — Xi(t)\ would always vanish in the 2r-periodic state. 

Next, we shall investigate conditions upon the coupling parameters C and r allowing for limit cycle oscillations. On 
one hand, if r becomes smaller than some T^m, the impulse from the excitatory neuron arrives too early to trigger a 
spike, since the system is still early in its refractory phase. On the other hand, if C becomes too small, the coupling 
force of the excitatory neuron is too weak to excite the system above its threshold and pull it far enough towards B. 

In Fig. [To] the regime of oscillations is shown in the parameter plane of the coupling strength C and coupling delay 
T. The oscillation period is color coded. The boundary of this colored region is given by the minimum coupling 
delay Tmm as a function of C. For large coupling strength, Tmin is almost independent of C; with decreasing C it 
sharply increases, and at some small minimum C no oscillations exist at all. At the boundary, the oscillation sets 





FIG. 10: Regime of oscillations in the (r, C) parameter plane for initial conditions corresponding to single-pulse-excitation in 
one system. The oscillation period T is color coded. The transition between black and color marks the bifurcation line. Inset 
(a) shows the oscillation period vs. r in a cut at C = 0.8. Inset (b): schematic plot of the saddle-node bifurcation of a stable 
(red solid line) and unstable (blue dashed) limit cycle. The maximal oscillation amplitude is plotted vs. the delay time r and 
the stable fixed point is plotted as a solid black line. The grey background marks the bistable region. Parameters: a — 1.05, 
e = 0.01. 
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in with finite frequency and amplitude as can be seen in the insets of Fig. [10] which show a cut of the parameter 
plane at C = 0.8. The oscillation period increases linearly with t. The mechanism that generates the oscillation is 
a saddle-node bifurcation of limit cycles (see inset (b) of Fig. [T0|) . creating a pair of a stable and an unstable limit 
cycle. The unstable limit cycle separates the two attractor basins of the stable limit cycle and the stable fixed point. 



5. DELAYED SELF-FEEDBACK AND DELAYED COUPLING 

In this section we consider the simultaneous action of delayed coupling and delayed self-feedback. Here we choose to 
apply the self-feedback term symmetrically to both activator equations, but other feedback schemes are also possible. 



eixi = Xi--^-yi + C[x2{t - r) - xi{t)] + K[xi(t - tk) - xi{t)] 
yi = xi+a 

€2X2 = X2-^ -y2+ C[xi{t - r) - X2{t)] + K[x2{t - Tk) - X2(t)] 

2/2 = X2+a (16) 

By a linear stability analysis similar to Sect. 4.1 it can be shown that the fixed point remains stable for all values 
of K and tk in case of a > 1, as without self- feedback. Redefining ^ = 1 — a^ — C — -fC, one obtains the factorized 
characteristic equation 

1 - + eA^ = XKe-^^'' ± ACe"^^ (17) 

Substituting the Hopf condition \ = iuj and separating into real and imaginary parts yields for the imaginary part 

-£_=K cos{ujtk) ± C cos(wt) (18) 

This equation has no solution for a > 1 since |^| ^a^ — 1 + C + K>C + K. 

The adopted form of control allows for the synchronization of the two cells not only for identical values of t and 
Tk-, but generates an intricate pattern of synchronization islands or stripes in the control parameter plane (Fig. Illj) 
corresponding to single-spike in-phase and antiphase oscillations with constant interspike intervals, see also Fig. fl^ a)- 
(d). Further, for adequately chosen parameter sets of coupling and self- feedback control, we observe effects such as 
bursting patterns Fig.fTWf) and oscillator death Fig.fTWcV In addition to these effects, there exists a control parameter 
regime in which the self- feedback has no effect on the oscillation periods (shaded yellow). 
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FIG. 11: Influence of delayed self-feedback upon coupled oscillations. The mean interspike interval (ISI) is color coded in the 
control parameter plane of the self- feedback gain K and delay tk- White areas mark regimes of irregular oscillations where 
the ISI variance becomes large (> 0.01). Time series corresponding to points (a)-(f) are shown in Fig. 1121 Other parameters: 
a = 1.3, e = 0.01, C = 0.5, r = 3. 
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FIG. 12: Different modes of oscillation corresponding to different self-feedback parameters K, t (red solid lines: activators 
Xi{t), green solid lines: inhibitors yi(t)). (a), (b): Antiphase oscillations for (a) K = 0.05, ta' = 3 (period T = 6) and (b) 
K = 0.5, TA' = 2 (T = 2); (c), (d): In-phase oscillations for (c) K = 0.5, tk = 3 (period T = 3) and (d) K = 0.5, ta- = 1.5 
(T — 1.5); (e): Oscillator death for K = 0.9, tk = 0.9; (f): Bursting pattern for K — 0.5, tk = 3.2. Other parameters: a — 1.3, 
£ = 0.01, C = 0.5, T = 3. 



Fig. [TT] shows the control parameter plane for coupling parameters of the uncontrolled system in the oscillatory 
regime {C — 0.5 and r = 3). We observere three principal regimes: (i) Control has no effect on the oscillation period 
(yellow), although the form of the stable limit cycle is slightly altered (Fig. [T^a)). (ii) Islands of in-phase and 
antiphase synchronization (color coded, see Fig. [12] (b)-(d)). (iii) Oscillator death (black) Fig. [T2l (e)). 

Fig- [131 shows the average phase synchronization time as a function of the coupling delay t and self- feedback delay 
Tk for fixed K = 0.5. The bright straight rays at rational tk/t indicate long intervals during which both subsystems 
remain synchronized. A particularly long average synchronization time is found if the two delay-times are equal. 
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FIG. 13: Average phase synchronization time (color coded) in the control parameter plane of coupling delay r and self-feedback 
delay tk- Other parameters: a = 1.3, e = 0.01, C — 0.5, K = 0.5. 



6. CONCLUSION 

Our analysis has focussed on a model of two coupled neurons, which may be viewed as a network motif for larger 
neural networks. We have shown that delayed feedback from other neurons or self- feedback from the same neuron 
can crucially affect the dynamics of coupled neurons. In case of noise-induced oscillations in instantaneously coupled 
neural systems, time-delayed self-feedback can enhance or suppress stochastic synchronization, depending upon the 
delay time. This offers promising perspectives with respect to potential therapies of pathological neural synchrony as 
occurring, for example, in Parkinson's disease. It suggests that by carefully choosing the delay time, feedback control 
applied locally to a neural sub-population can suppress the global synchronization of the neurons. 

In case of delay-coupled neurons without driving noise sources, the propagation delay of the spikes fed back from 
other neurons can induce periodic oscillations for sufficiently large coupling strength and delay times. Bistability of 
a fixed point and limit cycle oscillations occur even though the single excitable element displays only a stable fixed 
point. The two neurons oscillate with a phase lag of tt. If self-feedback is applied additionally, for example by axonal 
reflections (Fig. [T] (a)) in networks of electrically coupled pyramidal cells [IJI, synchronous zero-lag oscillations can 
be induced in some ranges of the control parameters, while in other regimes antiphase oscillations or oscillator death 
as well as more complex bursting patterns can be generated. 
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